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Abstract 

We establish the global existence and uniqueness of classical solutions to the 
Cauchy problem for the isentropic compressible Navier-Stokes equations in three 
spatial dimensions with smooth initial data which are of small energy but possibly 
large oscillations with constant state as far field which could be either vacuum or 
non-vacuum. The initial density is allowed to vanish and the spatial measure of 
the set of vacuum can be arbitrarily large, in particular, the initial density can 
even have compact support. These results generalize previous results on classical 
solutions for initial densities being strictly away from vacuum, and are the first for 
global classical solutions which may have large oscillations and can contain vacuum 
states. 



1 Introduction 



The time evolution of the density and the velocity of a general viscous isentropic 
compressible fluid occupying a domain 17 C is governed by the compressible Navier- 
Stokes equations: 

pt + div(pu) = 0, 

{pu)t + div(/9n ®u)- pAu - {p. + A)V(divu) + VP(p) = 0, 

where p > 0, u = {u^,u'^,v?) and P = ap^{a > 0,7 > 1) are the fluid density, 
velocity and pressure, respectively. The constant viscosity coefficients p and A satisfy 
the physical restrictions: 

/u>0, /i + ^A>0. (1.2) 
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University of Hong Kong. The research of J. Li is partially supported by NSFC Grant No. 10971215. 
Email: xdhuang@ustc.edu.cn (X. Huang), ajingli@gmail.com (J. Li), zpxin@ims.cuhk.edu.hk (Z. Xin). 



1 



Let = and p be a fixed nonnegative constant. We look for the solutions, 
{p{x , t) , u{x , t)) , to the Cauchy problem for (jl.ip with the far field behavior: 



u{x,t) — )■ 0, p{x,t) — 7- p > 0, as |x| — >• oo 



(1.3) 



and initial data, 



{p,u)\t=o = {po,uo) 



X e M^. 



(1.4) 



There are huge literatures on the large time existence and behavior of solutions to 
(|1.1|) . The one-dimensional problem has been studied extensively by many people, 
see |9ll2 H[3ni32j and the references therein. For the multi-dimensional case, the local 
existence and uniqueness of classical solutions are known in 128^133] in the absence 
of vacuum and recently, for strong solutions also, in [3l[5l[6l[3n] for the case that the 
initial density need not be positive and may vanish in open sets. The global classical 
solutions were first obtained by Matsumura-Nishida |27j for initial data close to a non- 
vacuum equilibrium in some Sobolev space . In particular, the theory requires that 
the solution has small oscillations from a uniform non- vacuum state so that the density 
is strictly away from the vacuum and the gradient of the density remains bounded 
uniformly in time. Later, Hoff [lOllllj studied the problem for discontinuous initial 
data. For the existence of solutions for arbitrary data (the far field density is vacuum, 
that is, /5 = 0), the major breakthrough is due to Lions [26] (see also Feireisl [7]), where 
he obtains global existence of weak solutions - defined as solutions with finite energy - 
when the exponent 7 is suitably large. The main restriction on initial data is that the 
initial energy is finite, so that the density vanishes at far fields, or even has compact 
support. However, little is known on the structure of such weak solutions. Recently, 
under the additional assumptions that the viscosity coefficients p and A satisfy 



and for the far field density away from vacuum (p > 0), Hoff ( |12tll4t [T5]) obtained 
a new type of global weak solutions with small energy, which have extra regularity 
information compared with those large weak ones constructed by Lions ( [26]) and 
Feireisl ( [7J). Note that here the weak solutions may contain vacuum though the 
spatial measure of the set of vacuum has to be small. Moreover, under some additional 
conditions which prevent the appearance of vacuum states in the data, Hoff ( |12y i5j) 
obtained also classical solutions. 

It should be noted that in the presence of vacuum, the global well-posedness of 
classical solutions and the regularity and uniqueness of those weak solutions ( [7 1ll2p26j ) 
remains completely open. Indeed, this is a subtle issue since, in general, one would not 
expect such general results due to Xin's blow-up results in [33], where it is shown that 
in the case that the initial density has compact support, any smooth solution to the 
Cauchy problem of the non-barotropic compressible Navier-Stokes system without heat 
conduction blows up in finite time for any space dimension, and the same holds for the 
isentropic case (jl.ip . at least in one-dimension, and the symmetric two-dimensional 
case ( [H]). See also the recent generalizations to the cases for the non-barotropic 
compressible Navier-Stokes system with heat conduction ( (5) and for non-compact 
but rapidly decreasing at far field initial densities ( [29]). 

In this paper, we will study the global existence and uniqueness of classical solutions 
to the Cauchy problem for the isentropic compressible Navier-Stokes equations, (jl.ip . in 
three-dimensional space with smooth initial data which are of small energy but possibly 



p > max{4A, —A} 



(1.5) 
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large oscillations with constant state as far field which could be either vacuum [p = 0) 
or non- vacuum > 0); in particular, the initial density is allowed to vanish, even has 
compact support. 

Before stating the main results, we explain the notations and conventions used 
throughout this paper. We denote 

fdx = / fdx. 
Jr3 

For 1 < r < cxD and /3 > 0, we denote the standard homogeneous and inhomogeneous 
Sobolev spaces as follows: 

'U = L'-(M3), zjfc-'- = {u€ Ll^iM.^) I \\V''u\\Lr < oo} , \\u\\Dk,r = WV^uUr, 
W^'"- = U r\ D^''' , H^ = W^^^, L)^ = L>'='2, = {ue L^\\\Vu\\l2 <oo} 



where / is the Fourier transform of /. 
The initial energy is defined as: 



Co = j Qpol^oP + G{po)^ dx, (1.6) 



where G denotes the potential energy density given by 



p 



It is clear that 



(G{p) = ^P, if p = 0, 

\ciip,p)ip-pf <Gip)<C2{p,p){p-p?, if p>0,0<p<p, 

for positive constants ci{p,p) and C2{p,p). 

Then the main results in this paper can be stated as follows: 

Theorem 1.1 Assume that \1.S\) holds. For given numbers M > (not necessarily 
small), f3 G (1/2,1], and p> p+1, suppose that the initial data {pq,uq) satisfy 

po\uo\^ + G{po) e L\ uoeH^nD^nD'', [p^- p,p{p^) - p{p)) ^ h\ (1.7) 

< inf po < suppo < P> WuidWijii^M, (1.8) 
and the compatibility condition 

- pAuo - (/U + X)Vdivuo + VP(po) = Po9, (1-9) 

for some g G with p^^g G L^. Then there exists a positive constant e depending on 
p, A, p, a, 7, p, (3 and M such that if 

Co < e, (1.10) 
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the Cauchy problem U.l]) il.^J]) ^-41 ) ^^ls a unique global classical solution (p, u) in 
X (0, oo) satisfying for any < r < T < oo, 

< p{x, t) <2p, xe M.^, t > 0, (1.11) 

'(p-p,P-P{p))GCi[0,T];H^), 
u G C{[0, T];D^n D^) D L'^{0, T; D^) n L°°(t, T; D^), 
ut £ L°°{0,T;D^)nL^{0,T;D^)nL°°{T,T;D^)nH^iT,T;D^), 

and the following large-time behavior: 

lim {\p- pi'' + p^/'^\u\^ + \Vu\^){x,t)dx = 0, (1.13) 

t— >oo J 

for all 

^ f(2,oo) /„.p>0, 
[(7,00), forp = 0. 

Similar to our previous studies on the Stokes approximation equations in |25j . we 
can obtain from (jl.lSp the following large time behavior of the gradient of the density 
when vacuum states appear initially and the far field density is away from vacuum, 
which is completely in contrast to the classical theory ( [151127] ). 

Theorem 1.2 In addition to the conditions of Theorem \l.l\ assume further that there 
exists some point xq G M'^ such that pq{xq) = 0. Then if p > 0, the unique global 
classical solution (p, u) to the Cauchy problem M.l]) 1^1.4^ obtained in Theorem 

li.il has to blow up as t —)• 00, in the sense that for any r > 3, 



lim ||Vp(-,t)||L. =00. 



A few remarks are in order: 



Remark 1.1 The solution obtained in Theorem ] becomes a classical one for positive 
time. Although it has small energy, yet whose oscillations could be arbitrarily large. In 
particular, both interior and far field vacuum states are allowed. 

Remark 1.2 In the case that the far field density is away from vacuum, i.e., p > 0, 
the conclusions in Theorem li.il generalize the classical theory of Matsumura-Nishida 
( ^0 c'^'SC of large oscillations since in this case, the requirement of small energy, 
il.lO\) . is equivalent to smallness of the mean-square norm of {po — p,uo). However, 
though the large-time asymptotic behavior /il.l3\) is similar to that in \27^ , yet our 
solution may contain vacuum states, whose appearance leads to the large time blowup 
behavior stated in Theorem li.^L this is in sharp contrast to that in ]15\ \2T^ where the 
gradients of the density are suitably small uniformly for all time. 

Remark 1.3 When the far field density is vacuum, i.e., p = 0, the small energy as- 
sumption, il.lO\} . is equivalent to that both the kinetic energy and the total pressure 
are suitably small. There is no requirement on the size of the set of vacuum states. In 
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particular, the initial density may have compact support. Thus, Theorem \1.1\ can be re- 
garded as a uniqueness and regularity theory of Lions- FeireisVs weak solutions in ^\26^ 
with small initial energy. It should also be noted that the conclusions in Theorem 
for the case of p = are somewhat surprising since for the isentropic compressible 
Navier-Stokes equations any non-trivial one- dimensional smooth solution with 

initial compact supported density blows up in finite time ( I34f ), and the same holds 
true for two-dimensional smooth spherically symmetric solutions ( 

Remark 1.4 It should be emphasized that in Theorem the viscosity coefficients 
are only assumed to satisfy the physical conditions While the theory on weak 

small energy solutions, developed in flE, [T5 ^ . requires the additional assumption lil.5\) 
which is crucial in establishing the time-independent upper bound for the density in the 
arguments in fT^. [T5 \l . 

Remark 1.5 For the incompressible Navier-Stokes system, a lot of results on the global 
wellposedness in scaling invariant spaces are available fR \22\. \23^I. In particular, Fujita- 
Kato 18] and Kato 122] proved that the system is globally wellposed for small initial 
data in the homogeneous Sobolev spaces H^^^ or in L^. In our case, the initial energy 
is small, therefore, we need the boundedness assumptions on the H^-norm of the initial 
velocity. It should be noted here that ^ L^/i^-m and 6/(3 - 2/3) > 3 for /? > 1/2, 
which implies that, compared with the results in f^.\22\l. our conditions on the initial 
velocity may be optimal under the smallness conditions on the initial energy. 

Remark 1.6 Similar ideas can be applied to study the case on bounded domain. This 
will be reported in a forthcoming paper 119] . 

We now comment on the analysis of this paper. Note that for initial data in the class 
satisfying (jl.7p - (|1.9p except uq G , the local existence and uniqueness of classical so- 
lutions to the Cauchy problem, (]l.ip - (|1.4p . have been established recently in [5]. Thus, 
to extend the classical solution globally in time, one needs global a priori estimates on 
smooth solutions to (ll.ip -( fL4]l in suitable higher norms. Some of the main new diffi- 
culties are due to the appearance of vacuum and that there are no other constraints on 
the viscosity coefficients beyond the physical conditions (II. 2p . It turns out that the key 
issue in this paper is to derive both the time-independent upper bound for the density 
and the time-depending higher norm estimates of the smooth solution {p,u). We start 
with the basic energy estimate and the initial layer analysis, and succeed in deriving 
an estimate on the spatial weig hted -norm of the velocity, the weighted spatial mean 
estimates on both the gradient and the material derivatives of the velocity. This is 
achieved by modifying the basic elegant estimates on the material derivatives of the ve- 
locity developed by Hoff ( [IO l ll2 U 13]) in the theory of small energy weak solutions with 
non- vacuum far fields and an interpolation argument. Then we are able to obtain the 
desired estimates on L^{0, min{l, T}; L°°(IR^))-norm and the time-independent ones on 
L8/3(niin{l,r},r; L~(R3)) -norm of the effective viscous flux (see ( 12.51) for the defini- 
tion). It follows from these key estimates and Zlotnik's inequality (see Lemma [2. 4p that 
the density admits a time-uniform upper bound which is the key for global estimates of 
classical solutions. This approach to estimate a uniform upper bound for the density 
is motivated by our previous analysis on the two-dimensional Stokes approximation 
equations in [25]. The next main step is to bound the gradients of the density and 
the velocity. Motivated by our recent studies ( 116 ^ 117 ^ 120]) on the blow-up criteria of 
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classical (or strong) solutions to such bounds can be obtained by solving a log- 

arithm Gronwall inequality based on a Beal-Kato-Majda type inequality (see Lemma 
I2.5P and the a priori estimates we have just derived, and moreover, such a derivation 
yields simultaneously also the bound for L^(0,T; L°°(R^))-norm of the gradient of the 
velocity, see Lemma 13.61 and its proof. It should be noted here that we do not require 
smallness of the gradient of the initial density which prevents the appearance of vacuum 
( |15tl27j). Finally, with these a priori estimates on the gradients of the density and 
the velocity at hand, one can estimate the higher order derivatives by using the same 
arguments as in [20] to obtain the desired results. 

The rest of the paper is organized as follows: In Section 2, we collect some elementary 
facts and inequalities which will be needed in later analysis. Section 3 is devoted to 
deriving the necessary a priori estimates on classical solutions which are needed to 
extend the local solution to all time. Then finally, the main results, Theorem 1.1 and 
Theorem 1.2, are proved in Section 4. 



2 Preliminaries 



In this section, we will recall some known facts and elementary inequalities which 
will be used frequently later. 

We start with the local existence and uniqueness of classical solutions when the initial 
density may not be positive and may vanish in an open set. 

Lemma 2.1 ( [5]) For p > 0, assume that the initial data (po > 0,uq) satisfy |i. 7p - 
l[1.9\) except uq G H^. Then there exist a small time T^: > and a unique classical 
solution {p, u) to the Cauchy problem (EiP ( lj.3|) (fi.^[ ) on R'^ x (0, T*] such that 

'ip-p,P-P{p))GC{[0,n];H'), 

uGCi[0,n];D^nD^)nL\0,n;D^), 
^ utGL^{0,n;D^)nL^{0,n;D^), e L°°(0, T,; L^), 

' ^UtteL^{0,T,;L''), t^/^uG L^{0,n;D^), 

t^l'^^putt G L°°(0,r,;L2), tut G L°°(0, T,; Z)3)^ 

tutt e L°°(o,r,;Z)i)nL2(o,r,;L>2). 

Next, the following well-known Gagliardo-Nirenberg inequality will be used later 
frequently (see [21]). 

Lemma 2.2 (Gagliardo-Nirenberg) Forp G [2,6],g G (l,oo), andr G (3, oo), there 
exists some generic constant C > which may depend on q,r such that for f G -ff^(M^) 
and g G ^^(M^) n D^'''{R^), we have 

ll/lli.<^^ll/ll£-^^/'llv/llir'^/^ (2.2) 



llnll / r'llnll'?(^-3)/(3^+'?(^-3))||V7„||3^/(3'-+<?(^-3)) 

\\9\\c{W) - l|V5llLr . [2.6) 

We now state some elementary estimates which follow from (j2.2p and the standard 
LP-estimate for the following elliptic system derived from the momentum equations in 

AF = div(/9u), pAuj = V X (pii), (2.4) 
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where 

f^ft + u-Vf, F^i2f, + X)divu-P{p) + P{p), oj^Vxu, (2.5) 
are the material derivative of /, the effective viscous flux and the vorticity respectively. 

Lemma 2.3 Let {p,u) be a smooth solution of il-l]) il.3\} . Then there exists a generic 
positive constant C depending only on p and A such that for any p £ [2,6] 

||VF||lp + \\Vuj\\lp < C\\pu\\lp, (2.6) 

WFWlp + II^IIl. < C\\pu\\fr'^^/^'^^ {\\Vuh^ + IIP - PmL^f-"^^^'"^ , (2.7) 
\\Vu\\lp <C(||F||iP + ||w||Lp) + C||P-P(^)||LP, (2.8) 
WVuUp < C||V.x||£-^)/('^) iWpuh- + IIP - P{p)he fP-'^/^'^^ . (2.9) 

Proof The standard LP-estimate for the elliptic system (j2.4p yields directly (j2.6p . 
which, together with (H^D and ([23]) . gives (f277l) . 

Note that —Au = — Vdivu + V x w, which implies that 

Vu = -V(-A)"^Vdivn + V(-A)"^V x co. 
Thus the standard estimate shows that 

llVnIliP < C(||divn||iP + ||i:i;||ip), for p G [2, 6], 



which, together with ([23]), gives Now ^M) follows from (I22D, and (|23]l . 

Next, the following Zlotnik inequality will be used to get the uniform (in time) upper 
bound of the density p. 

Lemma 2.4 ( [35]) Let the function y satisfy 

y'{t)=g{y) + h'{t) on [0,r], y{^) = y\ 
with g G C{R) and y,b e W^'^{0,T). If g{oo) = -oo and 

bit2) - Kti) < No + Ni{t2 - ti) (2.10) 
for all < ti < t2 < T with some Nq > and Ni > 0, then 

y{t) < max{2/°,C} +A^0 < oo on [0,r], 
where C, is a constant such that 

g{0<-Ni for C>C- (2.11) 

Finally, we state the following Beal-Kato-Majda type inequality which was proved 
in [1] when divu = and will be used later to estimate ||Vn||Loo and ||Vyo||/^2pii6. 

Lemma 2.5 For 3 < q < oo, there is a constant C{q) such that the following estimate 
holds for all Vu G L'^{m?) n D^^'i{M?), 

I|Vm||l-'(m3) < C (||divn||ioo(K3) + ||a;||ioo(iR3)) log(e + || V\||ig(K3)) 
+ C\\Vu\\l2(^z-.+C. 
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Proof. The proof is similar to that of (15) in [Ij and is sketched here for completeness. 
It follows from the Poisson's formula that 



\x - y\ 



dy 



= J divu{y)K{x - y)dy - j K{x - y) x Lo{y)dy 



V + w, 



where 



satisfies 



K{x - y) 



x-y 



A'k\x — yP ' 



\K{x -y)< C\x - yl-"^, \VK{x - y)\ < C\x - y\ 



(2.13) 



It suffices to estimate the term Vw since can be handled similarly (see [1]). Let 
6 E (0, 1] be a constant to be chosen and introduce a cut-off function ris{x) satisfying 
r]s{x) = 1 for |x| < 6,r]s{x) = for \x\ > 26, and |Vr/5(x)| < C6~^. Then Vv can be 
rewritten as 



= / 7]siy)K{y)Vdivu{x - y)dy - / Vr]s{x - y)K{x - y)divu{y)dy 



(2.14) 



+ y (1 - r]s{x - y))VK{x - y)dwu{y)dy. 
Each term on the righthand side of (j2.14p can be estimated by (j2.13p as follows: 

r]s{y)K{y)Vdivu{x - y)dy 



< C\\m{y)K{y)\\L,/(,-i)\\V'u\\L', 

/ /■25 \ (9-1)/? 

<C( r-2«/(«-i)r2drj 

< C7(5('?-=^)/'?||V\||l,, 



Vr]s{x - y)K{x - y)divu{y)dy 

< j \Vr,s{z)\\K{z)\dz\\divu\\L^ 

1-25 

<C J J^V^V^drlldivnllioo 

< C||divn||Loo , 

(1 - r]s{x - y))VK{x - y)d\\u{y)dy 



<C[ + \VK{x-y)\\d\wu{y)\dy 

\Js<\x-y\<l J\x-y\>\j 
-1 / /.oo \ 1/2 



(2.15) 



(2.16) 



(2.17) 



< Cy r-'^r^dr ||divn||Loo ^ C yj r-^r'dr 

< -Cln5||divti||Loo +C||Vn||i2. 



div?i||£,2 
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It follows from (^n^ - (^Tl\\ that 

\\Vv\\lo- < C (5'-1'^'>/''\\V'^u\\l9 + (1 - ln(5)||divn||Loc + \\Vu\\l2'^ . (2.18) 

Set 6 = mill 1 1, ||V2ti||^,'/('^"^)|. Then (fZTSD becomes 

||Vt>||L°° < C{q) (l + ln(e + ||V^n||L5)l|divn||Loo + ||Vu||i2) . 
Therefore ([232]) holds. 

3 A priori estimates 

In this section, we will establish some necessary a priori bounds for smooth solutions 
to the Cauchy problem (jl.ip (jl.3p ()1.4p to extend the local classical solution guaranteed 
by Lemma [2. II Thus, let T > be a fixed time and {p, u) be the smooth solution to (|l.ip 
(jl.SP (|1.4p on M'^ X (0,r] in the class (j2.ip with smooth initial data {po,uo) satisfying 
(|1.7p - (|1.9p . To estimate this solution, we set a{t) = min{l,t} and define 

Ai{T) = sup (a||Vti||^2) + / / ap\u\'^dxdt, (3.1) 
te[o,T] Jo J 



A2{T)= sup / p\u\'^dx+ / a^\Vu\'^dxdt, (3.2) 
te[o,T] J Jo J 

and 

AsiT) = sup / p\uf{x,t)dx. 

0<t<T J 

We have the following key a priori estimates on {p,u). 

Proposition 3.1 Under the conditions of Theorem \l.l{ for 

5o ^ (2/3 -l)/(4/3)e (0,1/4], (3.3) 

there exists some positive constant e depending on p, X, p, a, 7, p, 13 and M such that 
if {p, u) is a smooth solution of lll.3\) fil.4\ ) on x (0, T] satisfying 

sup p<2p, Ai{T) + A2{T) < 2C'J\ As{a{T)) < 20^^^^ , (3.4) 

M3 X [0,T] 

the following estimates hold 

sup p<7p/i, AiiT) + A^iT) < C'J\ AsiaiT)) < C','\ (3.5) 

M3x[0,T] 

provided Cq < £■ 

Proof. Proposition 13.11 is an easy consequence of the following Lemmas 13.31 13.41 and 
[331 

In the following, we will use the convention that C denotes a generic positive constant 
depending on p, A, p, a, 7, p, j3 and M, and we write C(a) to emphasize that C depends 
on a. 

We start with the following standard energy estimate for {p, u) and preliminary 
bounds for Vti and pii. 
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Lemma 3.1 Let {p,u) be a smooth solution of Ji.^P on x (0,T] with 

< p{x,t) < 2p. Then there is a positive constant C = C{p) such that 

sup [ (l-p\u\'^ + G{p)] dx + [ [ (p\Vu\'^ + {\ + p){dwuf)dxdt<Co, (3.6) 

0<t<T J V2 / JQ J 

Ai{T)<CCq + C j a\Vu\^dxdt, (3.7) 



/o 

and ^ 

A2{T) < CCo + CAi{T) + C j (T^\Vu\^dxdt. (3.8) 

Proof. Multiplying the first equation in (II. ip by G'{p) and the second by and 
integrating, applying the far field condition (jl.Sp , one shows easily the energy inequality 
M. 

The proof of ([SZD and I^M) is due to Hoff [lO]. For m > 0, multiplying l^^2 by 
C7™n and then integrating the resulting equality over lead to 

j a'^p\u\^dx = j (-a™n • VP + pa'^/^u ■u + {\ + ^)a'"Vdivu • u)dx 

3 (3.9) 

i=l 

Using (jl.ip i and integrating by parts give 



Ml = - / a™tt • VPdx 



= j (a"(divu)t(P - P(p)) - a'^iu ■ Vu) ■ VP)dx 

= (^j a"^divn(P - P{p))di!^ - ma'^-^a' j divn(P - P{p))dx 

+ y a™ (pV(divn)2 - P(divn)2 + Pd^u^dju^^ dx (^-lO) 

<(^y" cT"^divn(P - P(p))dx^ +mfT™'V'||P- P(p)||i2||Vn||i2 
+ C7(^)||Vtx||i2 

< a'"divn(P - P{p))dx^ + C{p)\\Vu\\l, + C(p)mV2('"-i)<T'C7o. 
Integration by parts implies 
M2 = y pa"^Au ■ iidx 

= [<y'"\\^Ah)t + ^^^'""'^^'llVnllia - /icT™ y ain^ai(n'=9fcn^)(ix (3.11) 



< _^ (cT^llVnll^^)^ + Cma™"^||Vn||22 + cy a"'\Vufdx, 
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and similarly, 

Ms = KlldivnlliO, + ^^^^^^--idivnlli. 

- (A + ^)cr" j divndiv(n • Vu)dx (3.12) 
< (CT^IIdiv-ull^a)^ + C7ma'"-i||Vn||2 2 + C J a'^lVufdx. 

Combining (13:91) - (I3l2]) leads to 

(cT™S(t))' + J (j'^p\u\^dx 

< {Cma"'-^ + C{p))\\Vu\\l2 + C {p)m^ a^'^"'-^'^ a' Cq + C j a'^\Vu\^dx, 



(3.13) 



where 



B{t) ^ ^W^uWl, + iA±if) lldivtxili. - j divn(P - P{p))dx 



> ^llV^lli. + ^^lldivnlli. - CC',/'\\di^uh. (3.14) 

>^\\Vn\\l, + ^-^\\divu\\l,-CCo. 

Integrating (|3.13|) over (0,T), choosing m = 1, and using (|3.14p . one gets (j3.7p . 

Next, for m > 0, operating a™"ii^ [d / dt + div(M-)] to (jl.ip n. summing with respect to 
J, and integrating the resulting equation over M.^, one obtains after integration by parts 

^ I p\u\'dx^ -!^a^-^a' I p\u\^dx 



(3.15) 



a^'ii^idjPt + div{djPu)]dx + p j a'^u^lAui + div(nAu.')]dx 
+ (A + /i) y (T'"'u.^ [5t5jdivu + div(ii9jdivn)]dx 



3 

A 



i=l 

It follows from integration by parts and using the equation (|l.ip i that 
iVi = - y a^'ii^idjPt + div{djPu)]dx 

= j a"'[-P'pdwudjU^ + dk{djU^u^)P-Pdj{dkU^u'')]dx (g jg) 

<C{p)a'^\\VuU2\\Vu\\L^ 
< Sa^'WVuWl^ + C{p,6)a"'\\Vu\\l2. 
Integration by parts leads to 

N2 = pj a^'ii^lAui + div{uAu^)]dx 

= -p I CT™[|Vttp + diii^dku'^diu^ - diii^diu'^dku^ - diU^ diu^ dku^]dx (3.17) 
a'^\Vu\'^dx + C j (T"|Vn|^dx. 



11 



Similarly, 

iV3<-^ 

6 - 2 

Substituting (j3.16p - ()3.18p into ()3.15p shows that for 5 suitably small, it holds that 



J a"'{divu)^dx + C J a'^lVul^dx. (3.18) 



(T^y p\u\'^dx^ +fij fT™|V'up(ix+ + A) J a"'{dwu)^dx 
< ma""^a' j p\u\^dx + Ca'^WVuWl, + C{p)a'^\\Vu\\l2. 



(3.19) 



Taking m = 3 in ()3.19p and noticing that 



we immediately obtain ()3.8p after integrating ()3.19p over (0,T). The proof of Lemma 
13. H is completed. 

Next, the following lemma will play important roles in the estimates on both Ai{(T{T)) 
[i = 1,3) and the uniform upper bound of the density for small time. 

Lemma 3.2 Let {p, u) he a smooth solution of ( li. ij) ^.3\) |i.^[ ) on x (0, T] satisfying 
Then there exist positive constants K and both depending only on p, \, p, a, 
7, p, j3 and M such that 

j-(T(T) j- 

sup t^-f^\\SIu\\l2+ t^-^ p\u\'^dxdt<K{p,M), (3.20) 

0<t<a{T) Jo J 

sup t'^-^ [ p\ufdx+ I ^ ^2-/3 I \\/u\^dxdt < K{p,M), (3.21) 

0<t<a{T) J Jo J 

provided Cq < Eq. 

Proof. As in [13], we define wi and W2 to be the solution to: 

Lwi = 0, wi{x,0) = wio{x), (3.22) 

and 

Lw2 = -VP(/5), W2{x,0) = 0, (3.23) 
respectively, with L being the linear differential operator defined by 

{Lwy = pwl + pu ■ Vw^ — {pAw^ + (p + X)divwxj) 
= pw^ - {pAw^ + (^ + X)divwx^), j = 1, 2, 3. 

Straightforward energy estimates show that: 

sup [ p\wifdx+ / ^ ^ /" \Vwi\'^dxdt < C{p) [ \wiofdx, (3.24) 

0<t<a{T) J Jo J J 

sup j p\w2\^dx+ j I \Vw2\^dxdt <C{p)Co. (3.25) 

0<t<a{T) J Jo J 
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It follows from (|3.22|) and standard L^-estimate for elliptic system that 

llVwill^e < C\\V'^wi\\l2 < C\\pwi\\l2. (3.26) 

Multiplying (j3.22p by wu and integrating the resulting equality over M^, we get by 
(|3:26D and ([33]) 3 that 

^{fi\\Vwi\\l2 + {fJ. + X)\\divwi\\l2)^ + j p\wi\'^dx 
pwi{u ■ Vwi)dx 

1/2 / /> \ 1/3 



<C{p)yj p\wi\^dxj [ I p\u\'^dx] ||Vwi||i6 
< C(p)CJ"/' / p\wi\''dx, 



which, together with Gronwall's inequality and (|3.24|) . gives 

,<t(T) 

livwiiii2+ / 

0<t<<T(T) 

and 

ra(T) 

0<t<o-(T) 



/>cr(T) I- 

sup ||Vu;i||^2+/ p\tdi\'^dxdt < CWVwioWli, (3.27) 

^t<<T(T) Jo J 

ra{T) I- 

sup t\\Vwi\\l2+ t p\wi\^dxdt<C\\wio\\l2, (3.28) 

:t<a(T) Jo J 



provided Cq < eoi = (2C(p))-3/'5o. 

Since the solution operator wio 1— )■ wi{-,t) is linear, by the standard Stein- Weiss 
interpolation argument ( [2J), one can deduce from (|3.27p and (j3.28p that for any 



l-cj{T) r- 

sup t^'^\\Vwi\\\2 + I t^'^ I p\wi\'^dxdt < C\\wiq\\\b, 

t<a(T) Jo J 



(3.29) 



0<t<a{T) 

with a uniform constant C independent of 6. 

Next, we estimate W2- It follows from a similar way to ()2.6p and (|2.8p that 

f \\Vii2p + A)divt/;2 -{P- Pimh^ < C\\pw2\\l^, 

\\\Vw2\\Le<Ci\\pW2\\L^ + \\P-Pip)\\L^)- 

Multiplying (j3.23p by W2t, integrating the resultant equation over and using (j3.30p . 
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one has 
1 
2 



^(^fl\\Vw2\\l2+{fl + X)\\dwW2\\l2-2 j {P - P{p))divw2dx^ + j p\td2\^dx 
J pw^iu • Vw2)dx — j Ptdivw2dx 

C{p) (^j p\w2\^dx^ ' y p\ufdx^ ' \Vw2\\l^ 

+ y divt(;2div((P - P{p))u)dx + j {P{p) + (7 - l)P)divndivt(;2(ix 



< C{p)cl''^ (1 p\w2\^dxy (\\p^'^W2\\l^ + \\P-PmL^\ 
- j{P- Pip))n ■ V (div«;2 - dx 

+ 2{2p + X) J^^~ + C\\Vu\\l2 + C\\Vw2\\l2 

< C{-p)cl''^ j p\w2\''dx + CCy"" + C\\P - P{p)U^ \\uU4p"^W2\\l^ 
+ C\\P-P{p)\\\, + C\\Vu\\l2+C\\Vw2\\l. 

<C(p)Co^°/' / p\w2?dx + CCl'^ + C\\Vu\\l2 + C\\Vw2\\l2, 



which, together with p.25p and Gronwall's inequality, gives 



sup \\Vw2\\l2 + / / p\w2?dxdt < CCl^^, (3.31) 

t<a(T) Jo J 



0<t<a(T) 

provided Co < £02 = {2C{p))-^^^<-\ Taking 

■u^io = ^^0 so that wi + W2 = u, we then 
conclude from (fH^OD and (fH^TT) that for any 6 £ [/3, 1], 

rcr{T) I- 

sup t^"lVn||22+/ t^~^ p\u\'^dxdt<C\\uQ\\\,+CCl'^, (3.32) 

0<t<tj{T) Jo J 

provided Co < eo — niinjeoi, £02}- Thus, p.20p fohows from p.32p directly. 

To prove (I3.2ip . we take m = 2 — /? in (I3.19P to obtain, after integrating (I3.19P over 
(0,(7(T)) and using (l3:32]l and dZl]), that 

sup t''^ [ p\u\^dx+ r^^^2-/3 f 
0<t<(7(T) J Jo J 



7n| dxdt 



l-cryi ] 

<C t^-^WVuW^^idt + Cip.M) 

Jo 



<C I t'~^\\Vu\\L2 {\\pu\\l2 + IIP - P{p)\\h) dt + C{p, M) 
'^{T) , ,1/2 

^(2/3-l)/2 (tl-/3||Vn||i.) ' (t2-/3||pV2^||2^)l/2^^1-/3||^l/2^||2^)^^ 



< c 

+ C{-p,M) 



<C(p,M)( sup t^^^ /"p|n|2(ix ) +C(p,M) 

\o<t<o-(T) J y 
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which imphes (j3.2ip . Thus, we finish the proof of Lemma |3.2[ 
The following Lemma 13.31 will give an estimate on A3(a{T)). 

Lemma 3.3 // {p, u) is a smooth solution of /il.3\) on x (0, T] satisfying 

Ii3.4^ , there exists a positive constant ei depending on p, A, p, a, 7, p, j3 and M such 
that the following estimate holds for defined by i3.3\) : 



sup 

0<t<a{T) 



p\u\^{x,t)dx <Cq°, (3.33) 



provided Cq < ei. 



Proof. Multiplying (jl.ip ? by 3|n|u, and integrating the resulting equation over M^, 
we obtain by ()2.9p that 

Jt J ^'"'''^^ 

<C j \u\\Vu\'^dx + C j \P - P{p)\\u\\Vu\dx 

< C\\u\\l4Vu\\%^\\Vu\\]^,^ + C\\P - P{p)\\l3\\u\\l6\\Vu\\l2 

< C\\Vuff i\\pu\\L2 + IIP - PiMi^f^ + cc'J''\\Vu\\l, 

1/2 



< C'WVuff (\\pu\\l^ + + CCo/'||Vn||| 



lL2 



<Ct(2'5o-3/2)(l-/3)(il-/3||y^||2^)-25o+5/4(^l-/3||^l/2^||2^)l/4||y 

+ CCo'/^V3(i~'^)/4(ti-/3||Vtx||2,)3/4||Vn||i2 + CC',/''\\Vu\\l„ 
which together with (I3.20|) and ()3.6p gives 

sup / p\u\'^dx 
0<t<o-(T) J 

<C{p,M)[ / t 3-8*0 dt] / \\Vu\\l2dt] 



|4(5o 
Il2 



, / MT) \ / \ 1/2 (3.34) 

+ C{-p,M)Cl'''{j^ t-3(i-/^)/2dt n ||Vn||i.dt 

+ j po\uQf'dx + CCQ 
<C{-p,M)Cl'\ 

provided Cq < eo; where in the last inequality we have used the following simple facts: 

/ r \ 3(2/3-l)/(4/3) 

po\uo\^dx <C PQ\uQ\'^dxj ll^^oll^^f'^^ 
<C{-p,M)Cl'\ 

and 

2(3 -4Jo) (1-/3) _ /3(2/3-l) 
3 - 8(5o 2-13 
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due to (1331) and /3 G (1/2, 1]. Thus, it follows from (pm|l that (|CTD holds provided 
Co < El, where 

ei ^ min {eq, (C(p, M))-i/'5o | = ^^^^ ^c{p, M))-^^/(2/3-i) | _ 
The proof of Lemma 13.31 is completed. 

Lemma 3.4 There exists a positive constant e2{f^, p, a, 'j, p, P, M) < Si such that, if 
{p, u) is a smooth solution of il.3\) ^1-4^ on x (0, T] satisfying then 

AiiT) + A2{T) < Cl^\ (3.36) 

provided Cq < £2- 

Proof. Lemma |3 . 1 1 shows that 

Ai{T)+A2{T)<C{p)Co + C{p)[ a^\\Vu\\l,ds + C{p) [ aWVuWl^ds. (3.37) 

Jo Jo 

Due to (USD, 

[ a^\\Vu\\l4ds<C [ {\\F\\% + WojWI^) ds + C [ a^\\P - P{p)\\l4ds. (3.38) 
Jo Jo Jo 

It follows from ([221) that 

r a'iWFWl, + Ml) ds 
Jo 

<C a' {\\Vuh2 + IIP - P{p)\\l^) WpuWhds 
Jo 

< C(p) sup fa3/2||^n||^2 (ai/2||Vn||i2 +Co/^)) / [ ap\u\'^dxds 
te{o,T] ^ ^ ^ ^ Jo J 



< C{p) [A^iT) + Co^/^j ^^/^(T)yli(T) 

< C(p)Co. 

To estimate the second term on the right hand side of (I3.38P , one deduces from (jl.ip 1 
that P — P{p) satisfies 

(P - P{p))t + u ■ V(P - P{p)) + 7(P - P(^))divtt + 7P(^)divn = 0. (3.40) 

Multiplying (j3.40p by 3(P — P{p))'^ and integrating the resulting equality over M'^, one 
gets after using divu = ^^^^ {F + P — P{p)) that 

= -[j{P- P{~P)fdx^^ - 0^ j{P- P{p)fFdx 

-3jP{p) j {P - P{p)fd\wudx 
<-l^j{P- P{p)fdx^+ r^\\P - P{p)\\\, + C7^||F||i4 + C,||Vn||i.. (3.41) 
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Multiplying (|3.41|) by a^, integrating the resulting inequality over (0, T), and choosing 
rj suitably small, one may arrive at 



a'\\P-Pmhdt 



< C sup IIP - Pip)\\h +C IIP - Pkhjul^ 

0<t<T Jo 



a(T) 

^ dt 



(3.42) 



+ C{p) r a^\\F\\l,ds + C{p)Co 
Jo 

< C{p)Co, 

where p.39p has been used. Therefore, collecting ()3.38p . ()3.39p and ()3.42p shows that 

l^a^ {\\Vu\\% + IIP - P(p)||i4) ds < Cip)Co. (3.43) 
Jo 

Finally, we estimate the last term on the right hand side of (I3.37p . First, (I3.43P 
implies that 

/ / a\Vufdxds < [ [ (|Vn|^ + \Vu\'^)dxds < CCq. (3.44) 

Ja{T) J Ja{T) J 

Next, one deduces from ([23]), (13:20]) and ([33]) that 
a\\Vu\\'l3dt 

,-a(T) 

ra(T) ^ ^ ^ iff X^''^ 

^11/2 ( t I p\u\^dx j ( 



< C{p) £ t\\Vu\\%^ {\\pu\\l2 + Cq/') dt 

<C{p)J [t^^-^^/^\\Vu\\L2j\\Vu\\]^2i^tJp\u\'^dxj dt + C{p)Co 
<C{p) sup U'^-^^''^\\Vu\\l2) r^^\\Vu\\y^ (t f p\u\'^dx\ 

iG(0,(7(T)l ^ ^ Jo \ J J 



(3.45) 



E(0,<7(T)] 

+ C(p)Co 
<C{-p,M)A^^^Cl'* + C{p)C, 
<C{-p,M)Cl'\ 

provided Cq < £1. It thus follows from (|3.37p and (|3.43p - (|3.45p that the left hand side 
of ()3.36p is bounded by 

C{-p,M)Cl'' <Cl'^ 

provided 

Co<e2 = min{ei,(C(^,M))-s}. 

The proof of Lemma 13.41 is completed. 

We now proceed to derive a uniform (in time) upper bound for the density, which 
turns out to be the key to obtain all the higher order estimates and thus to extend the 
classical solution globally. We will use an approach motivated by our previous study 
on the two-dimensional Stokes approximation equations ( [25j). 
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Lemma 3.5 There exists a positive constant e = e{p,M) as described in Theorem \l.l\ 
such that, if (p, u) is a smooth solution of lil.l]) on M'^ x (0, T] satisfying 

iT4\ ), then 

7 — 

sup < -T' 

0<t<T 4 

provided Cq < e. 

Proof. Rewrite the equation of the mass conservation as 

Dtp = g{p)+h'{t), 

where 

For t e [0,ct(T)], one deduces from Lemma [221 UM^ dSSS]), (fCTj) . ([MT]) and flO]) 
that for 5o as in ()3.3p and for ah < ti < t2 < '^{T), 

\b{t2)-h{t^)\ 

<C \\{pF){;t)\\L^dt 



<C{p) r \\F{.M!'\NF{-M\!edt 



< C{-p) / \\p"^u\\^l^\\Vu\\^lUt 
Jo 



< c{p) £^''\-('-^y^\\pu\\];,' (t'-f'wvuwi.y^' dt 

[ t-^'-^y'\\pu\\fdt\ 

lrr\ \ 3/4 



dt 



= C{p,M) t-K2-m-So+2/3)+5o] (^^2-/3||^l/2^||2^^-^°+l/3 (^t\\p'/2^2^^^ 

<C{-p,M){A^{a{T))f<^l^ 
< C{p,M)cf'''^, 

provided Cq < £2- Therefore, for t G [0,o"(T)], one can choose Nq and A''i in (|2.10p as 
follows: 

AT, = 0, No = C{p,M)Cl^''\ 
and C = /5 in (fZTT]) . Then 

9{0 = -lz^{Q^ -p')<-Ni = Q, for all C > C = /5. 
Lemma 12.41 thus yields that 

sup IIpIIl- < max{^, p] + < p + C{p, M)cf < ^, (3.46) 

te[o,(7(T)] ^ 
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provided 

Co < min{e2,e3}, fores 



8/(3<5o) / - N 32/3/(3(2/3-1)) 



2C{p,M)J \2C{p,M) 



On the other hand, for t G [a{T),T], one deduces from Lemma [221 (I3.36P . (I3.6p . and 
2M that for all a{T) <ti<t2<T, 

\h{t2)-h{h)\<C{p) r \\F{;t)\\Lo^dt 
Jti 



< :r^it2 -h) + C{p) [ \\F{.,t)\fi!'dt 

<^-{t2-h) + C{-p)Cl''' r \\Vui;t)\\hdt 



2/i + A 



provided Cq < £2- Therefore, one can choose A''i and Nq in ()2.10p as: 



Note that 



So one can set C = P + 1 in (|2.1ip . Lemma 12.41 and p.46p thus yield that 

sup WpWl^ < max I ^, /5 + i| + iVo < V + C{p)C^o^'' < ^, (3.47) 

provided 

Co <e = m.m{e2, 63,64}, for £4 = i ^^^-s^ j ■ (3-48) 

The combination of (j3.46p with ()3.47p completes the proof of Lemma 13.51 

From now on, we will always assume that the initial energy Cq satisfies ()3.48p and 
the positive constant C may depend on 

1 /2 

T, \\Po 9\\l^, ||Vc/||i2, ||Vno||//2, \\po-p\\h3, \\P(po) - P{p)\\h3, 

besides p, A, p, a, 7, p, (3 and M, where 5 is as in (|1.9p . 

Next, we will derive important estimates on the spatial gradient of the smooth solu- 
tion (p, u). 

Lemma 3.6 The following estimates hold 



sup / p\u\'^dx+ [ [ \Vu\'^dxdt <C, (3.49) 

0<t<T J Jo J 

sup (||Vp||L2nL6 + IIVuIIhO + / \\Vu\\L--dt<C. (3.50) 

<t<T Jo 
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Proof. Taking ^ = 1 in (|3.32|) together with (j3.36p gives 

te[o; 



sup ||Vii||^2 + / / p\u\'^dxdt < C. (3.51) 
;efo,Ti Jo J 



Taking m = in (I3.19p . one can deduce from Gaghardo-Nirenberg's inequality (12. 2p . 
(USD, (|33T1) and (fXT9]) that 

j p\u\'^dx^ +p j |Vn|^(ix + (/i + A) j {diwiifdx 



<C\\Vu\\\, + C{-p)\\Vu\\l, 
< C\\Vu\\l2\\Vu\\Ic.+C 



L6 

<C{\\F\\le + \\u\\l, + \\P-Pmh)+C 

< C(||VF|||2 + \\Vuj\\l,)+C 

< C\\pu\\l2+C 
<C\\p^/Mi2+C. 

Taking into account on the compatibility condition (ll.9p . we can define 

^u{x, t = 0) = ^/p^g. (3.53) 

Then (j3.49p follows from (|3.5ip - (j3.53p and Gronwall's inequality. 

Next, we prove (j3.50p by using Lemma 12.51 as in [17] . For 2 < p < 6, |Vp|^ satisfies 

{\Vp\P)t + div(|V/)|fn) + {p- l)|V/)|fdivn 

+ p\yp\p''^{ypYyu{vp) + pp\vp\p~'^vp ■ vdi^u = o. 



(3.52) 



Thus, 



due to 



dtWVpUv < C(l + ||Vn||ioo)||Vp||LP + C\\V^u\\lp 
< C(l + ||Vn||Loo)||V/>||LP + C\\pu\\lp, 



(3.54) 



||V\||lp < C{\\pu\\lp + ||VP||lp) , (3.55) 
which follows from the standard L^-estimate for the following elliptic system: 

pAu + {p + A)Vdivti = pu + VP, u — )• as |x| — )• oo. (3.56) 

It follows from Lemma 12.51 and p.55p that 

||Vu||lcx) < C (lldivnilioo + ||a;||Lcx)) log(e + ||V^ii||j^6) + C||Vii||j;^2 + C 

< C (lldivnllioo + ||a;||L=o) log(e + \\u\\i6 + ||VP||^6) + C 

< C (||divn||ioo + llwllicx)) log(e + ||V'u||/^2) 
+ C (lldiviillLoo + ||w||ioo) log(e + || VpH^e) + C. 

Set 

fit) ^ e + ||Vp||i6, g{t) ^ 1 + (lldivnlli^ + ||l^||l^) log(e + \\Vu\\l2) + \\Vu\\l2. 
Combining ()3.57p with ()3.54p and setting p = 6 in ()3.54p , one gets 
fit) < Cg{t)f{t) + Cg{t)f{t)\nf{t) + Cg{t), 



(3.57) 
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which yields 

{Infit))' <Cg{t) + Cg{t) In fit), (3.58) 
due to /(t) > 1. Note that ([23]), Lemma [221 <^M: (1339]), and Lemma [33] imply 

rT fT 

/ g{t)dt<C {\\divu\\lc^ + \\uj\\loo) dt + C 
Jo Jo 

<c [ (||F||ioo + IIP - p(/5)||ico + ||w||ioo) dt + c 

Jo 

<C I {\\F\\l2 + \\VF\\le + IMl^ + llVwIlie) + C 
Jo 

< C / ||Vn||^2dt + C 
Jo 

<C, 

which, together with ()3.58p and Gronwall's inequality, shows that 



(3.59) 



Consequently, 



sup fit) < C. 

0<t<T 



sup WVpWlh < C. (3.60) 

0<t<T 



As a consequence of ()3.57p . ()3.59p and (j3.60p . one obtains 

rT 



\Vu\\Lo^dt < C. (3.61) 







Next, taking p = 2 in (j3.54p . one gets by using (j3.6ip . (j3.5ip and Gronwall's inequality 
that 

sup \\Vp\\l2 < C, 

0<t<T 

which, together with P3SD . ([232]), (l?3T]l . <^M), and gives ([530]) . The proof 

of Lemma [3.61 is completed. 

The following Lemmas 13.7113.10] will deal with the higher order estimates of the so- 
lutions which are needed to guarantee the extension of local classical solution to be a 
global one. The proofs are similar to the ones in [20|, and we sketch them here for 
completeness. 



Lemma 3.7 The following estimates hold 

rT 

0<t<T J Jo 



sup ! p\ut\^dx+ ! [ \Vut\^dxdt <C, (3.62) 

<t<T J Jo J 



sup {\\p - p\\h2 + \\P{p) - P(p)||^2) < C. (3.63) 

tG[0,T] 
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Proof. Estimate (|3.62p follows directly from the following simple facts: 



< C + C\\p^/^u\\l2\\u\\l6\\Vu\\16 



and 



|V«t|||2 < llVttll^a + ||V(n • Vu)\\l2 



< ||Vtt||^2 + ClliillioollV^nll^a + C||V'u||24 

< \\Vu\\l2 + C, 

due to Lemma [3161 

Next, we prove (I3.63p . Note that P satisfies 

Pt + u-VP + 'jPdivu = 0, (3.64) 

which, together with (jl.ip i and a simple computation, yields that 

IdlV'Plli. + IIVVlli.) (3,3, 

< C7(l + ||Vn||ioo) (IIV^PII^^ + llvVllia) + C||F||^2 + CMl^ + C, 

where we have used the following simple fact: 

||Vu||/fm < C (||diVti||j7m + ||w||//m) 

<C{\\F\\Hm + \\uj\\Hm + \\P - P{p)\\h-^), for m=l,2. 



Noticing that F and uj satisfy (I2.4p . we get by the standard L^-estimate for elliptic 
system, ([OO]) and ([330]) that 

11-^11^2 + ||w||h2 < C'(II^IIl2 + \\u\\l2 + \\pu\\l^ + \\V{pu)\\l2) 
<C{1 + \\Vp\\l3\\u\\l^ + \\Vu\\l2) 
< C(l + ||Vn|U2), 

which, together with (j3.65p . Lemma 13.61 and Gronwall's inequality, gives directly 

sup {\\V^P\\l2 + ||VV||l2) < c. 
te[o,T] 

Thus the proof of Lemma 13.71 is completed. 
Lemma 3.8 The following estimates hold: 

sup i\\pt\\m + \\Pt\\m)+ [ {\\pa\\l2 + \\Ptt\\l2)dt<C, (3.66) 

0<i<T Jo 



sup j\Vut\^dx+ f [puldxdt<C. (3.67) 

0<t<T J Jo J 
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Proof. We first prove (j3.66p . One deduces from (j3.64p and (|3.5U|) that 

WPth^ < C||n||Loo||VP||i2 +C||Vn||i2 < C. 
Differentiating p.64p yields 

VPt + u ■ VVP + Vu ■ VP + 7VPdivu + 7PVdivn = 0. 
Hence, by ()3.50p and (|3.63p . one gets 

||VPt||i2 < C\\u\\lo^\\V'^P\\l2 +C||Vn||i3||VP||i6 +C\\V'^u\\l2 < C. 
The combination of ()3.68p with p.69p imphes 

sup \\Pt\\m < C. 

0<t<T 

Note that Pu satisfies 

Ptt + -fPtdivu + -fPdivut + ut ■ VP + u • VPt = 0. 
Thus, one gets from ([3771]) ([3770]) ([330]) and iK62\i that 



(3.68) 



(3.69) 
(3.70) 

(3.71) 



rT 



dt 



<C [ (||Pt||L6||V^z||i3 + \\Vut\\L2 + WutheWVPhs + WVPth^fdt 

Jo 

< C. 

One can handle pt and pu similarly. Thus (I3.66P is proved. 

Next, we prove ()3.67p . Differentiating (jl.ip ? with respect to t, then multiplying the 
resulting equation by uu, one gets after integration by parts that 



■ ^ ^-^ jpt\ut\'^dx - jptu ■ Vu ■ utdx + jPtdiYUtdx 
+ ^ jpu\ut\'^dx + j{ptu ■ Vu)t ■ Utdx - jpUfVu- uudx 



(3.72) 



pu ■ Vut ■ Uudx — I Pttdivutdx 

5 



i=l 



It follows from ^Jifi, (ICTIl . I^Mh and (IHTU^D that 



pt\ut\'^dx — ptu-Vu-utdx+ Ptdwutdx 



< 



1 
2 

div {pu)\ut\'^dx 



+ C'll/OtllLsllu • Vu||i2||Ui||^6 + C||Pt||i2||VUt||i2 



<C yp|u||ut||Vut|dx + C||Vnt||i2 

< C||7x||i6||pi/2^t||J^/'||nt||^/'||Vnt||i2 + C||V^xt||i2 
<<5||Vnt|||2 + Q, 



(3.73) 
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2|/i 



Ptt\ut\ dx 



{ptu + put) ■ V{\ut\'^)dx 



and 



<C\\Vnt\\l,+C\\Vut\\%' 
<C\\Vut\\l2 + C, 

J [ptu ■ Vu)^ ■ utdx 

jipuu.Vu.ut + ptUt.Vu.u, + p,u.Vut.ut)dx 

< IIPmIIlHI'" • Vullisllutll/^e + ||/5t|lL2|ll^tl^llL3||Vn||i6 
+ ll/^tllL^II^IlL"" II Vut 11^2 11^6 

<C||p«||i2+C||V^Zi||i2. 



L2 



Cauchy's inequality gives 



I/3I + I/4 



put • Vu ■ uttdx 



+ 



pu ■ Vut ■ Uttdx 



< C\\p^^'^utt\\L2 (||nt||i6||Vti||i3 + ||n||Loo||Vnt||i2) 

< dWp'/^uttWh + CsWVutW 



L2' 



and 



Pttdivutdx 



< ||Ptj||i2||divnt||i2 
<C\\Ptt\\l. + \\Vutf 



L2- 



Due to the regularity of the local solution, (j2.ip . tVut G C([0, T*]; L^). Thus 

||V^/t(-,r,/2)||i2 < |-||tVtXt||ioo(o,T.;L2) 
J- * 



(3.74) 



(3.75) 



(3.76) 



(3.77) 



(3.78) 



where C may also depend on ||V5||j^2. 

Collecting all the estimates one deduces from I^7l7\i . (fSrB^j) and 

Gronwall's inequality that 



T,/2<t<T 

On the other hand, ()2.ip gives 



sup \\Vut\\L2 + [ [p\ 

\/2<t<T Jn/2J 



Uttl'^dxdt < C. 



(3.79) 



sup ||V'Uj||i2 

0<t<T,/2 



T,/2 



p\utt\^dxdt < C. 



(3.80) 



The combination of (|3.79p with (|3.80p yields (|3.67p immediately. This completes the 
proof of Lemma [ 
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Lemma 3.9 It holds that 



sup (IIp-pII^s + ||P-P(^)||^3) < C, (ogi) 
te[o,T] ^ ' ' 



sup (||V'Ui||i2 + ||Vu||^2)+ / (||V'u||2^3 + ||V'Ut||2^i)dt<C. (3.82) 

SfO.Tl Jo 



fT 

(||V'Ui||i2 + ||VU||^2) + / 

te[o,T] 

Proof. It follows from (l3:671l and (1330)1 that 



(3.84) 



||V(H||i2 < \\\Vp\\ut\\\L2 + ||pVnt||i2 + |||V/9||^z||Vtx|||i2 

+ \\p\^A^\\l^ + Ilpkliv'^x|||i2 

< ||Vp||i3||nt||i6 +C7||Vnt||i2 +C7||Vp||i3||n||L^||Vn||i6 
+ C|| Vn||/,3 ||Vti||^6 + C||u||lcx) ||V^ti||^2 

<c, 

which together with ()3.49p gives 

sup ||pn||/^i < C. (3.83) 

0<i<T 

The standard ff^-estimate for elliptic system (j3.56p gives 

||V^n||j:^i < CII^An + (/i + A)Vdivti||j:^i 
= C||pn + VP||j:^i 
<C(||/5^||^i + ||VP||^i) 
<C, 

due to (|l.ip 9. (j3.83p and (|3.63p . As a consequence of (j3.50p and (|3.84p . one has 

sup ||Vu||j:^2 < C. (3.85) 

0<i<T 

Therefore, the standard L^-estimate for elliptic system, (|3.5Up . and Lemma ElSl yield 
that 

||V^ut||/^2 < C||/iAnt + (^ + A)Vdivnt||i2 

= \\puu + Ptut + Ptu ■ Vn + put ■Vu + pu- Vut + VPtlUa 

< C{\\pUu\\l^ + ||/0t||L3|kt||L6 + ll/0t||L3|l^i|U°°l|V'"llL6) (3.86) 
+ C{\\ut\\L4^u\\L3 + \\u\\L^\\Vut\\L2 + ||VPt||i2) 

< C\\pUtt\\L2 + C, 

which, together with (j3.67p . implies 

rT 

llVntll^idt < C. (3.87) 

Applying the standard ff^-estimate for elliptic system (|3.56p again leads to 

WV'^uWh^ < C\\pAu + {p + A)Vdivn||j|^2 

< C\\pu\\h^ + C||VP||h2 (3.88) 

<C + C\\Vut\\Hi+C\\V^P\\L^, 



25 



where one has used (j3.83p and the fohowing simple facts: 

\\^Hput)\\L^ < C (|||VVl|nt||U2 + |||VH|Vnt|||i2 + || WH^^) 

< C [\\V^p\\L2\\Vut\\m + l|Vp||L3||Vtxt|U6 + II WIU2) 
<C + C\\Vut\\m, 

and 

\\V\pu-Vu)\\l2 < C(|||V2(p^x)||Vn|||i2 + |||V(HI|V'^x|||i2 + ||V\||i2) 

< + ||V2(pn)||L2||V^z||H2 + ||V(pu)||i3||V2n||i6) 

< C7 (1 + ||VVllL2hl|L- + l|V/9||i6||Vn||i3 + llV^nllia) 

due to (I3:63]l and ([335]). By usmg (IMS]), ([5:88]) . and (IMS]), one may get that 
(llV^PlliO, 

< C(|||V\||VP|||i2 + |||v2n||V2p|||i2 + |||Vn||V^P|||i2 + ||V^n||i2) ||V^P||i2 

< C(||V\||i2||VP||j^2 + ||V\||i3||V2p||i6 + ||Vu||Lco||V3p||i2) ||V^P||i2 

+ C{1+ \\V^Ut\\L2 + ||V3P||^2) ||V3p||i2 

<C + C\\Vut\\]j,+C\\V''P\\l,, 
which, together with Gronwall's inequahty and ()3.87p . yields that 

sup ||V^P||i2 < C. (3.89) 

0<t<T 



Collecting all these estimates (j3.87p - (j3.89p and (j3.63p shows 

0<t<T Jo 

It is easy to check similar arguments work for p — phy using ()3.90p . Hence, 



sup ||P-P(p)||j^3+ / \\Vufj^sdt<C. (3.90) 



sup ll/f) — pII j|^3 < C. (3.91) 

0<t<T 

Combing ([3:90]) with ^m\i shows l[3M\i . Estimate ([3:82]) thus follows from (13:67]) . 
([3:85]) . ([3:87]) . and ^Mh . Hence the proof of Lemma EJ] is finished. 

Lemma 3.10 For any r G (0,T), there exists some positive constant C{t) such that 
sup (||Vnt||j^i + ||V^u||i2) + / j\Vutt?dxdt < C{t). (3.92) 

T<t<T Jt J 

Proof. Differentiate (|l.ip 9 with respect to t twice to get 
puttt + pu ■ ^utt - pAutt - {p + A)Vdivntt 

= 2d\\{pu)utt + d\N{pu)tUt - 2{pu)t ■ Vut - {pttu + 2ptUt) ■ Vu (3.93) 
- putt ■ Vti - VPtt- 
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Multiplying (j3.93p by Uu and then integrating the resulting equation over M^, one gets 
after integration by parts that 



1 d 



jp\utt\'^dx + j (/i|Vnttp + + A)(divutt)^) dx 

4 julfpu ■ Vul^dx - j{pu)t ■ [V{ut ■ Uu) + 2Vut ■ utt] dx 

jipttu + 2ptut) ■ Vn • uttdx - jpuu ■ Vn • uudx + jPttdivuudx 



(3.94) 



5 

A 



1=1 

We estimate each Ji(i = 1, • • • , 5) as follows: 
Holder's inequality gives 

l-Zil < CWp^^'^uuWl^WVuuWl^WuWl^ 

<6\\Vuu\\h+Cs\\p'/^uu\\h. 

It follows from (^Mh . ([MZD, and ([330]) that 

\J2\ < C{\\pUt\\L3 + \\ptu\\L3) {\\uu\\L4^Ut\\L^ + ||VUtt||i2||nt||i6) 

< c(||/)^/V|li('hi|lii' + llptllL^hlUe) WVutth^ (3.96) 
<6\\Vutt\\l2+Cs, 



(3.95) 



l-^al < C {\\ptt\\L2\\u\\L^\\Vu\\L3 + \\pt\\L<i\\ut\\Le\\^u\\L'2) WuuWl^ 



< d\\Vuu\\l2 + CsWptt'^ 



L2 ^ <^S\\Ptt\\L2, 



and 



(3.97) 
(3.98) 



+ l^sl < C'llmt||L2||V'"llL3||^tti|lL6 + CWPuWl'^W^UuWl^ 

< 6\\Vuu\\l2 + CsWp'/^uttWh + CsWPttWh. 

For any r E (0,r*), since t^^'^^uu G L~(0,r*;L2) by ^Jii, there exists some io e 
(r/2, r) such that 

p\uu\'^dx{to) < ^||i^/^^/p««|li^(o,T,;L2) ^3 gg^ 

< C{t). 

Substituting (j3.95p - (j3.98|) into (|3.94|) and choosing 6 suitably small, one obtains by 
using ()3.66p p.99p and Gronwall's inequality that 



sup p\utt\^dx+ / hVuttl'^dxdt < C{t), 
to<t<T J J to J 

which, together with (j3.86p and (j3.67p . yields that 

sup WVutWm + r huu\^dxdt < C(r), (3.100) 

T<t<T Jt J 

due to to < T. Now, (l3:92]l follows from (l3:88]l . (I31MI1 . and (ISlHTTl . We finish the proof 
of Lemma 13.101 
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4 Proof of Theorems 11.11 and 11.21 



With all the a priori estimates in Section [3] at hand, we are ready to prove the main 
results of this paper in this section. 

Proof of Theorem By Lemma l2.ll there exists a T^, > such that the Cauchy 
problem (jl.ip . (jl.3p . (jl.4p has a unique classical solution {p,u) on M'^ x (0, T*]. We will 
use the a priori estimates, Proposition 13.11 and Lemmas 13.91 and 13.1^ to extend the 
local classical solution (p, u) to all time. 

First, it follows from ([3351) and ([L8]) that 

Ai (0) + A2 (0) = 0, ^3 (0) < Co^° , po<p, 

due to Co < e. Therefore, there exists a Ti G (0,T^,] such that (|3.4p holds for T = Ti. 
Next, we set 

T* = sup{r I (m holds}. (4.1) 

Then T* > Ti > 0. Hence, for any < t < T < T* with T finite, it follows from 
Lemmas 13.91 and 13.101 that 

Vut,V^u£C{[T,T];L^nL'^), Vu, V^u G C (^[r, T]; n C (^S^)) , (4.2) 

where we have used the standard embedding 

L°^{T,T;H^)nH\T,T;H-^) ^ C {[t,T]; L'^) , for any q G [2,6). 
Due to (IM2|), (IMZI), and (^Mi), one can get 




<c, 

which yields 

pi/VeC([r,T];L2). 
This, together with (j4.2p . gives 

p^/2n,VnG C7([r,r];L2). (4.3) 



Next, we claim that 

T* = oo. (4.4) 

Otherwise, T* < oo. Then by Proposition EH ([33]) holds for T = T*. It follows 
from Lemmas 13.91 and [3.101 and (14. 3p that {p{x,T*),u{x,T*)) satisfies (ll.7p - (ll.9p except 
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u{-,T*) G H^^ , where g{x) = u{x,T*), x G M'^. Thus, Lemma 12.11 imphes that there 
exists some T** > T* , such that ([O]) holds for T = T** , which contradicts ([0]) . 
Hence, ()4.4p holds. Lemmas 12. H 13.9] and 13.10) and (14.20 thus show that (p, u) is in fact 
the unique classical solution defined on M"^ x (0, T] for any < T < T* = oo. 

Finally, to finish the proof of Theorem II. H it remains to prove (jl.lSp . 
Multiplying ()3.40p by 4(P — P{p))^ and integrating the resulting equality over M^, 
one has 

{\\p - p{p)\\uy (t) 

= -(47 - 1) jiP- P{p))^divudx - j P{p){P - P{p)fd\-vudx, 



which yields that 

/OO 
(IIP-P 

due to ([333]). Combining (|3:43]l with (j43]) leads to 

lim IIP -P(/5) 11^4 =0 

r— >oo 

which together with ()3.6p implies 



+ ||V'u||^4) dt < C, 



lim / \p — pl'^dx = 0, 

t— i>00 J 



for ah q satisfying (fTTil) . Note that ([3:6]) and (H^]) imply 

1/2 



y p^/^l^dx < (^J p\u\'^dx^ \\u\\% < C\\Vu\\l2. 



Thus ()1.13p follows provided that 



lim ||Vti||j;^2 = 0. 

t— >oo 



Setting 



T/,\ A IIV7 ii2 , + I 
Ht) = ^\Nu\\l^ + I 



divnl 



L2, 



choosing m = in (|3.9p . and using (j3.1ip and (|3.12p . one has 

|/'(t)| <C j p\u\'^dx + CWVuWls + CCI^^\\Vu\\l2, 
where one has used the following simple estimate: 



|Mi| 



■ii • VPdx 



(P - P{p))divudx 



< CCq^^\\Vu\\l2. 



(4.5) 



(4.6) 



(4.7) 
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We thus deduce from (gTT]), ([233), and (IHIiSD that 

<C J {\\p^^^u\\l2 + ||Vn||^4 + llVttll^a) dt 

which, together with 

\I{t)\'^dt <C J \\Vu\\l2dt < C, 

imphes ()4.6p . The proof of Theorem II .11 is finished. 

Proof of Theorem M.SX The proof is similar to that of Theorem 1.2 in [25]. We just 
sketch it here. 

Otherwise, there exist some constant Ci > and a subsequence {tnj^^_^^ t^- 

OO such that || Vp (•, ||^^ < Ci. Hence, the Gaghardo-Nirenberg inequahty (|2.3p 
yields that there exists some positive constant C independent of t„. such that for 
a = r/(2r-3) G (0,1), 

\\p{x,tn^) - p\\c{W) 

< C ||V/)(X,^.)||2, \\p{x,tn^) - pW^'" 

<CCt\\p{x,tn^)-p\\^-\ (4.8) 
Due to (|1.13|) . the right hand side of (|4.8|) goes to as — ?■ oo. Hence, 

||p(x,t„J - /o||(^(-p-) as OO. (4.9) 

On the other hand, since (/>, n) is a classical solution satisfying (|1.12|) . there exists a 
unique particle path xo(t) with xo(0) = xq such that 

p{xQ{t),t) = for ah t > 0. 

So, we conclude from this identity that 

||p(2;,in,) -p||t7(K3) > \p{xo{tn^),tn^) - p\ = P > ^, 

which contradicts (|4.9|) . This completes the proof of Theorem 1.2. 
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